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1.  Introduction 


The  U.S.  Army's  soldiers  will  make  unprecedented  use  of  advanced  systems  and  sensors  on  the 
Future  Combat  System's  (FCS's)  battlefield  (see  http://www.darpa.mil/fcs/index.html, 

03/26/02).  This  will  include  the  use  of  low-cost,  non-l I ne-of-sight  acoustic  sensor  systems  for  the 
surveillance,  detection,  identification,  classification,  and  tracking  of  military  targets  (i.e.,  Fong 
and  Srour  1994;  Rosenthal  and  Stevens  1994;  Srour  and  Robertson  1995;  Y  oung  et  al.  1999; 
Becker  and  Gudesen  2000;  and  M  ays  and  Price  2000).  M  ost  new  U.S.  Army  acoustic  systems 
use  unattended  microphone  sensors  to  construct  small  ground-based  beam-forming  arrays  to 
determine  the  I ine-of- bearing  angle  of  sound-emitting  targets  [e.g.,  the  remote  netted 
acoustic/seismic  sensor  array  (RNADS)  (Thompson  and  Scherer  1991;  Carnes  and  M  organ 
1995;  Wilson  1997;  Eom  etal.  1999;  Lopez  etal.  1999;  Wellmann  and  Srour  1999;  Wellman 
1999;  Depireux  and  Shamma  2000;  and  M  ays  and  Vu  2000)].  Flence,  the  U.S.  Army  has  a 
growing  interest  to  implement  accurate  and  reliable  computer  models  for  determining  point-to- 
point  acoustic  transmission  (West  etal.  1991;  West  etal.  1992;  Wilson  1993;  Nobleand  Marlin 
1995;  and  Wilson  2000).  At  the  same  time,  the  retrieval  and  interpretation  of  acoustic  signals  in 
diverse  microclimate  area,  (e.g.,  in  and  around  forests,  hilly  terrain,  arid  deserts,  coastal  areas,  or 
in  cities)  is  greatly  influenced  by  turbulence  and  refraction  effects  caused  by  finer  scale 
atmospheric  motions  over  varying  topography  and  surface  energy  budgets  (Smith  1989; 
Auvermann  and  Goedecke  1993;  Auvermann  etal.  1995;  Goedecke etal.  1997;  Wilson  1998ab; 
Goedecke  etal.  2000;  and  Auvermann  2001).  We  expect,  therefore,  that  improved  physics- based 
theory  and  computer  models  for  meteorology  coupled  to  acoustics  will  contribute  important 
information  on  the  performance  of  critical  battlefield  systems. 


2.  Meteorology  and  Acoustics 


The  topic  of  this  report  is  the  coupling  of  meteorology  to  acoustics  in  forests.  We  focus  on 
deriving  estimates  for  sound  speed,  from  which  it  is  possible  to  determine  sound  intensity  and 
attenuation.  The  value  for  the  speed  of  sound  may  be  computed  from  Newton's  2nd  Law  as 
follows  (Fleagleand  Businger  1963):  The  change  in  pressure  {dp)  over  a  unit  cross-sectional  area 
of  a  fluid  containing  a  sound  wave  (i.e.,  between  the  plane  of  compression  and  plane  of 
rarefaction)  can  be  expressed  as 

dp  =  ds  —  ,  (1) 

dt 

where  is  the  density  of  the  fluid  medium,  is  an  increment  of  distance  in  the  direction  of 

propagation,  and  dc  is  the  change  in  sound  speed  over  an  increment  of  time,  dt.  Equation  (1)  can 
be  re-written  as 


-  dp  =  c  dc 


(2) 


1 


because  = 


-1 


ds 

is  specific  volume  and  —  =  c.  The  volume  (y)  of  fluid  moving  through  a  unit 

dt 

cross-sectional  area  of  a  fluid  medium  \sV  =  cdtso  that  any  fractional  change  in  this  volume 
(^)  can  be  equated  to  a  fractional  change  in  sound  speed  (i.e.,  ^  =  — ).  Also,  the  volume  of 

fluid  is  proportional  to  its  specific  volume  (i.e.,  V  =  M  ,  and  jy  =  M  J  ,)  which  yields 

—  =  —  =  —  .  Therefore,  equation  (2)  can  be  re-written  as 
y  c 


dp 


(3) 


Based  on  the  principles  of  thermodynamics  and  the  behavior  of  ideal  gases,  one  expects  that  this 
ratio  of  pressure  change  to  volume  change  will  depend  on  the  amount  of  heattransWed 
between  adjacent  fluid  elements  undergoing  either  compression  or  expansion.  Textbooks  such  as 
FI  eagle  and  Businger  (1963)  consider  the  foil  owing  two  extreme  cases:  Heat  transfer  is 
sufficiently  rapid  such  that  no  temperature  grad  lent  exists  in  the  fluid  (i.e.,  temperature  remains 
constant  during  the  process);  or  heat  flux  is  negligible  and  the  process  is  adiabatic.  To 

RT 

demonstrate  the  first  case  we  can  expand  the  equation  of  state  (EOS)  (i.e.,  p  = — )  where  R  = 

M 

8314.32  J  moF^  is  the  universal  gas  constant  and  M  is  molecular  mass,  for  constant 
temperature  (r)  as 

p  d  +  dp  =  0  ,  (4) 


SO  that  substitution  into  equation  (3)  yields, 

2  RT 

c  =  p  =  —  or  c 
M 


(5) 


This  sound  speed  for  isothermal  (constant  temperature)  expansion  is  often  called  the 
'Newtonian'  velocity  of  sound. 

In  the  second  case,  where  pressure  and  volume  change  across  elements  of  a  fluid  medium, 
transmitting  sound  waves  are  considered  to  occur  under  adiabatic  conditions.  Then  the  value  for 
the  speed  of  sound  may  be  computed  from  the  1st  Law  of  Thermodynamics  as  follows  (Fleagle 
and  Businger  1963): 

dH  =  dQ  +  p  d  ,  (6) 


where  JiT  is  the  incremental  amount  of  heat  added  to  a  given  system,  dQ  =  CvTis  the  internal 
energy  of  the  system,  where  Cv  is  the  specific  heat  at  constant  volume  of  a  fluid  (in  units  /  Kg~^ 
fC^),  and p  d  is  the  work  or  change  in  energy  of  the  system  as  it  undergoes  expansion  or 

compression.  Because  an  adiabatic  process  by  definition  is  one  in  which  no  heat  is  added  or  lost 
to  a  given  system  (i.e.,  dH=  0)  and  by  definition  an  isentropic  process  is  that  where  no  entropy 

dH 

is  added  or  lost  (i.e.,  dS^  =  —  =  0)  where  is  the  differential  of  specific  entropy,  then  from 
the  EOS  and  equation  (6)  we  can  derive  the  following  expression: 


2 


(7) 


0 


M 


Because  —  =  Cp  c^,  where  Cp  is  the  specific  heat  at  constant  pressure  of  a  fluid  (in  units  7^^ 
then  through  some  algebraic  manipulation,  equation  (7)  can  be  re-written  as, 


0  = 


Cp  d 
Cy 


(8) 


Returning  to  the  problem  to  define  the  ratio  of  pressure  change  to  volume  change  for  an  adiabatic 
(or  isentropic)  process,  we  have  from  equation  (8), 


dp  _  Cp  p 
d  Cy 


(9) 


so  that  substitution  into  equation  (3)  yields. 


c  = 


RT 


M 


(10) 


where  =  ^  is  the  ratio  of  specific  heats.  Equation  (10)  for  the  speed  of  sound  (c)  has  been 


found  to  agree  closely  with  observations.  For  application  to  outdoor  acoustics,  Wong  and 

Embleton  (1984, 1985)  have  deduced  the  ratio  of  specific  heats  and  molar  mass  ( — )  as  a 

M 

function  of  temperature  and  humidity  in  the  form. 


—  =  0.04833 +  (/z  0.023)A  ,  (11) 

M 

where  A,  =  9.2  10"^ -1-5.5  10"®? -i- 4.25  10"V, /z  is  humidity  (in  the  range  0.0  to  1.0)  and  ris 
air  temperature  in  degrees  Celsius.  From  equation  (11),  Wong  (1986)  has  calculated  the  value  for 
the  speed  of  sound  in  standard  dry  air  at  0  °C  and  a  barometric  pressure  of  1013.25  mbar  to  be 
Co  =  331.29  ms"^ 

In  addition,  to  account  for  increases  and  decreases  in  sound  speed  due  to  variations  in  wind 
velocity,  it  is  useful  to  define  an  effective  sound  speed  from  the  following  expression  given  by 
Noble  and  Marlin  (1995)  and  Osteshev  (1997): 

Cejf  =  -I-  M  COS  (  ^  s  p)  ,  (12) 
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where  M  isthemeanof  the  horizontal  wind,  is  the  bearing  of  the  wind  from  North,  /^isthe 

bearing  of  the  receiver  from  the  source,  and  u  cos  [  ^ is  the  component  of  the 

sound  speed  along  the  direction  of  propagation  from  the  source  to  the  receiver  (see  Figure  1).  As 
discussed  in  Osteshev  (1997),  the  effective  sound  speed  in  equation  (12)  is  valid  only  for  nearly 
horizontal  propagation  angles  (Wilson  2002). 


Figure  1.  Geometry  schematic  for  the  expression  Cgy  =  co  +  M  cos  (  -  s.r), 

where  Cejfls  the  effective  sound  speed  that  includes  the  sound  speed  component 
due  to  the  mean  wind  along  the  path  of  propagation. 


The  effective  speed  of  sound  in  equation  (12)  will  also  vary  with  height  (z)  above  ground  level  as 
a  function  of  the  profiles  of  air  temperature,  humidity,  and  wind  velocity.  Generally,  sound  speed 
differences  across  vertical  layers  will  cause  acoustic  waves  to  be  refracted  upward  if  the  effective 
sound  speed  decreases  with  height  and  refracted  downward  if  sound  speed  increases  with  height. 
Over  flat  and  uniform  terrain,  wind  velocity  profiles  through  the  surface  layer  (z<  50  m)  are 
generally  increasing  (logarithmically)  with  height  and  are  easily  predicted  via  traditional  profile 
theory  (M  onin  and  Obukhov  1954).  Likewise,  temperature  profiles  can  be  extrapolated  upward 
through  the  surface  layer  for  homogeneous  and  steady  state  conditions.  I  n  contrast,  computer 
models  to  describe  meteorological  profiles  within  and  above  forests  are  not  as  easily  determined 
or  well  known.  Therefore,  the  remainder  of  this  report  will  present  an  overview  of  selected  past 
research  on  meteorology  in  forests  to  examine  the  calculation  of  the  speed  of  sound  through  the 
atmosphere  in  the  forest  environment  for  military  acoustic  applications.  Our  objective  is  to 
evaluate  meteorological  models  for  estimating  wind  speed  and  temperature  profiles  within  and 
above  forests  to  determine  which  models  are  most  applicable  in  representing  mechanical  and 
thermodynamic  influences  on  the  speed  of  sound  in  the  forest  environment. 
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3.  Vertical  Profiles  of  Wind  Speed  in  Forests 


An  extensive  survey  of  literature  shows  that  the  study  of  winds  in  forests  has  been  motivated 
primarily  by  scientific  interests  to  better  understand  surface  exchanges  of  heat,  moisture  (water 
vapor),  and  carbon  dioxide.  Such  properties  not  only  influence  regional  and  local  climate,  but 
also  affect  the  forest  ecosystem,  growth,  and  regeneration.  Others  have  retrieved  wind  related 
data  in  forests  for  the  study  of  transport,  diffusion,  and  deposition  of  air-borne  pollutants  and 
aerial  sprays.  At  the  same  time,  scientists  have  found  it  necessary  to  assess  the  impact  of  wind 
damage  on  trees  and  tree  harvests  by  severe  storms  (e.g.,  Coutts  and  Grace  1995). 

3.1  Meteorological  Measurements 

Numerous  vertical  wind  speed  profiles  within  and  above  forest  canopies  have  been  retrieved 
through  meteorological  measurements  and  reported  by  various  authors  (e.g.,  Allen  1968;  Bergen, 
1971;  M artin  1971;  Oliver  1971;  Raynor  1971;  Denmead  and  Bradley  1985;  Amiro  and  Davis 
1988;  Baldocchi  and  Hutchinson  1988;  Baldocchi  and  M  eyers  1988;  and  Lee  and  Black  1995). 
The  main  features  of  wind  speed  profile  structure  within  and  above  a  forest  stand  (Figure  2)  on 
clear  days  are  as  follows:  1)  a  low-level  wind  maximum  at  0.1  to  0.3A,  where  the  live  canopy 
(i.e.,  branches  and  leaves)  extends  from  about  0.3h<z<h,  where  h  is  the  height  of  the  canopy 
top;  2)  a  layer  of  minimum  wind  speed  at  0.6  to  0.8A,  which  coincides  with  the  region  of 
maximum  leaf  area  density;  and  3)  a  diurnal  time  maximum  of  wind  speed,  which  extends 
through  the  entire  canopy  in  the  early  afternoon.  M  ini  mum  wind  speeds  through  the  tree  crown 
have  been  associated  with  the  zero-plane  displacement  height,  d.  Zero-plane  displacement 
determined  from  various  experiments  in  pine  forests  have  fallen  within  the  range  0.70  -  0.76A 
(Bergen  1971;  Raynor  1971;  Oliver  1971;  and  Lee  and  Black  1993).  Bergen  (1971)  remarked 
that  similar  features  in  profile  structure  have  resulted  from  airflow  measurements  through 
modeled  forests  in  wind  tunnels  (e.g.,  M  eroney  1968;  Sadeh  etal.,  1971;  and  Chen  et  al.,  1995). 
Oliver  (1971)  commented  that  over  the  course  of  his  fieldwork  no  significant  variations  in  profile 
shape  were  found  for  different  mean  wind  speeds. 
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Figure  2.  Typical  vertical  wind  speed  profile  structure  within  and  above  a 
forest  stand  (based  on  Bergen  1971). 


However,  in  a  recent  measurement  study  of  the  vertical  distribution  of  wind  speed  in  a  dense 
maize  canopy  {h  =  1.7m),  Jacobs,  etal.  (1995)  found  a  (sub-canopy)  local  velocity  maxima 
around  0.1/?.  Similarly,  Lee  and  Black  (1993)  reported  secondary  maxima  in  wind  speed  at  0.12/? 
fora  relatively  uniform  expanse  of  tall  pine  trees  (/?  =  16.7m,  on  average).  Hence,  the  height  of 
the  subcanopy  wind  maxima  is  strongly  influenced  by  the  shape,  density,  and  leaf  area 
distribution  of  the  trees. 


M  assman  (1982)  studied  vertical  distributions  of  foliage  (i.e.,  needle  surface  area,  for  nine  old- 
growth  Douglas  fir  crowns  and  one  old-growth  sugar  pine)  and  showed  several  comparisons  of 
data  and  distribution  models  to  estimate  the  observed  profiles.  Analysis  of  M  assman's  data  and 
other  leaf  area  models  reported  in  the  literature,  such  as  M  eyers  (1987),  Paw  U  and  M  eyers 
(1989)  and  M  eyers  etal.  (1998),  suggested  that  forest  canopies  regularly  conform  to  one  of  three 
general  leaf  area  distribution  profiles  (Figure  3).  They  found  that  leaf  area  distributions  are  not 
always  symmetric  about  the  layer  of  maximum  foliage  density,  like  Profile  1,  but  are  more  often 
skewed  upward  toward  the  top  of  the  forest  canopy,  as  shown  in  Profile  2  and  Profile  3.  In  turn, 
M  eyers  et  al.  (1998)  presented  useful  parameter! zations  for  roughness  length  (zo)  and  a  zero- 
plane  displacement  height  (J)  as  a  function  of  leaf  area  density,  based  their  previous  works  as 
well  as  the  earlier  work  of  Shaw  and  Periera  (1982), 


d  =  h 


0.05  + 


LAI 


0.2 


■  + 


jP  1) 
20 


(13) 


and 
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^0 


h  0.23 


(14) 


Lai^25 

10  67 

where  LA/ is  the  leaf  area  index  and  L  =  1,  2,  or  3  correspond  to  one  of  the  three  general  profile 
shapes  described  in  Figure  1. 


Figure  3.  Normalized  vertical  profiles  of  leaf  area  distribution  for  forest  canopies  where 
a(z)=  {z)  lai  (based  on  M  eyers  etal.  1998;  and  Albertson  etal.  2001). 


h 

By  definition,  the  leaf  area  index  isLA/=  A(z)  dz,  where  Afzj  is  the  leaf  area  density  through 

0 

the  small  vertical  layer  between  z  and  z  +  dz  per  unit  surface  area  of  ground  below  (M  unn  1966). 
V alues  for  leaf  area  index  for  forests  vary,  but  have  been  reported  most  often  in  the  range  LA/=1 
to  5  (Pielke  1984;  Kaimal  and  Finnigan  1994;  Finnigan  and  Brunet  1995;  and  Pyles  etal.  2000). 
Roughness  length  (zg) ,  in  contrast,  is  the  quantity  that  varies  in  the  following  manner:  zo  is 
-lO""*  m  over  snow,  sand,  dry  lakebeds,  or  concrete,  ~10"^m  over  soils,  ~0.1  m  over  farmlands, 
tall  grass,  or  shrubs,  and  ~1.0  m  for  orchards  and  forests  (Oke  1978  and  Pielke  1984).  Typical 
values  for  these  types  of  data  also  can  be  found  in  a  table  compiled  by  RacheleandTunick 
(1994). 

3.2  Wind  Speeds  Above  an  Idealized  (Uniform)  Forest  Stand 

Such  profile  characterizations  assume  (idealized)  steady  and  horizontally  homogeneous 
conditions  over  large  areas.  Initially,  therefore,  we  may  assume  that  at  some  height  above  tall 
roughness  elements,  wind  speed  profiles  are  similar  in  form  to  those  descri bed  for  lower  roughness 
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elements  (Busch  1973).  Asa  result,  wind  speed  profile  structure  above  an  Idealized  forest  stand 
can  be  described  as 


d  U  _  m 
d  z  Kz  d) 


(15) 


where  U  Is  the  mean  wind  speed,  m  =  Is  the  friction  velocity,  —  =  u'w'  Is  the  vertical 

turbulent  shearing  stress.  Is  air  density,  k  Is  Karman's  constant  {k  ~  0.4),  d  Is  zero-plane 
displacement,  and  ^  is  a  non-dimensional  parameter  that  is  a  function  of  atmospheric  (thermal 
and  mechanical)  stability  (M  onin  and  Obukhov  1954  and  Dyer  1974).  Then  based  on  the  work  of 
Panofsky  (1963)  and  Benoit  (1977),  the  mean  vertical  profiles  of  wind  speed  above  forests  can 
be  determined  as 


U  =  In  ^  ^  ,  (16) 

k  ZQ 

where  m  is  referred  to  as  thediabatic  influence  function,  which  modifies  the  standard  log-law 
wind  profile  for  non-adiabatic  lapse  rate  conditions  (Paulson  1970  and  Webb  1970). 

Alternately,  M  artin  (1971)  applied  a  power-law  exponent  model  to  determine  the  influence  of 
stability  and  site  characteristics  on  profile  structure  in  the  layer  above  the  canopy  top, 

,  p 

Uz  =Uz  .  (17) 

where  and  are  the  wind  speeds  (in  units  m/s)  at  heights  zi  and  Z2  above  the  canopy  top 

(A),  and  theexponentO  <P<  1  isdetermined  from  measured  wind  speed  data  for  different  times 
of  day  (i.e.,  for  varying  thermal  stability).  Several  authors  have  provided  tables  to  give  examples 
of  predicted  and  observed  values  for  the  power-law  exponent,  P,  for  various  sites  (DeM  arrais 
1959;  Touma  1977;  and  Irwin  1979). 

Immediately  above  the  forest  stand,  however,  there  Isa  region  within  which  modified  log-law 
and  power-law  exponent  models  tend  to  depart  from  the  observed  profiles  (Garratt  1992).  The 
so-called  transition  layer,  or  roughness  sub-layer,  transfers  additional  momentum  from  the  wind 
flow  towards  the  surface  due  to  canopy-top  generated  turbulence  (Garratt  1992;  Kaimal  and 
Finnigan,  1994;  and  Zoumakis  1994).  Garratt  (1980, 1992)  contends  that  atmospheric 
observations  from  numerous  field  studies  supports  the  following  empirical  modification  to  the 
wind  speed  and  wind  shear  profiles  for  tall  crops  and  trees 

3  JJ  ^  / 1  Q  \ 

“  77  A  772  m  • 

d  z  k(z  d) 

where 

*=exp  0.7  1  —  for  z<z*  ,  (19) 
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where  z  >  h  is  height  aboveground  level,  and  z*  is  the  depth  of  the  roughness  sub-layer.  Values 

of  —  for  stable  and  unstable  conditions  are  reported  to  vary  between  10  and  about  150  and  have 
zo 

been  found  to  be  a  function  of  roughness  density,  (i.e.,  lower  values  corresponding  to  more  dense 
canopies  and  vise  versa  [Garratt  1980, 1992]).  Asa  result,  the  depth  of  the  transition  layer  for 
typical  forests  (zo  ~1.0  m)  was  found  to  vary  according  to  z*  ~3  ,  where  is  tree  spacing 

(  =  10-30  m,  possibly).  A  Iternately,  A  rya  (2001)  suggests  that  the  top  of  the  transition  sublayer 

could  be  estimated  asz*  =  1.5A  to  2.5/?,  also  depending  on  the  roughness  density.  Equation  (18) 
is  integrated  to  yield  the  mean  vertical  profiles  of  wind  speed  through  the  transition  sublayer  as 

U  Uj^=  —  In  -  m  +  m  foi'  z  <  z*  ,  (20) 

k  zo 

where  =  /|  *  is  defined  in  the  same  manner  as  „,(i.e., 

*1  *1  2 

=  -i-ln  — ^  2tan  ^.x-i-y,  where  ;c  =  ^  see  (Paulson  1970)  and  (Webb 

1970).  ^  isa  positive  quantity  and  gives  the  excess  in  wind  speed  due  to  enhanced  canopy-top 

turbulence,  which  in  effect,  decreases  the  vertical  profile  gradient  (see  Figure  4).  In  a  calculation 
described  by  Zoumakis  (1994),  the  modeled  forest  height,  displacement,  and  roughness  were 
given  ash  =  18.5  m,  J  =  11.3  m,  and  zo  =  1.2  m,  respectively,  which  resulted  in  a  value  for  the 
depth  of  the  transition  layer  of  z*  =  32  m. 


Figure  4.  M  ean  vertical  profiles  of  wind  speed  through  the 
transition  sublayer  for  neutral  (adiabatic) 
conditions  using  equation  (20)  to  account  for 
enhanced  canopy-top  generated  turbulence.  The 
log-law  profile  is  modified  by  m,  which  gives 
an  excess  in  wind  speed  and  decreases  the 
profile  gradient.  For  this  example,  h  =  10  m,  zo 
~  0.74  m,  u*  =  0.5  ms‘\  and  =  30  m. 
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3.3  Wind  Speeds  Within  an  Idealized  (Uniform)  Forest  Stand 

From  conservation  of  momentum  for  steady  state,  homogeneous  (non-advective  and  equilibrium) 
conditions,  airflow  within  the  canopy  and  subcanopy  layers  of  an  idealized  (uniform)  forest 

stand  is  often  described  as  the  balance  of  momentum  flux  divergence  (^-^^-'^),  horizontal 

d  z 

pressure  gradient  (^-^),  and  drag  due  to  friction  (-QA  u^),  which  is  the  drag  imparted  to  the 
d  X  2 

wind  flow  by  the  canopy  leaves  and  branches  (Barr  1971;  Shaw  1977;  and  Shinn  1971).  The 
equation  for  wind  flow  within  the  forest  canopy,  therefore,  is  often  expressed  as 


d  u'  w' 
d  z 


d  p 

d  X 


\ciA 


(21) 


where  Q  = 


is  the  surface  drag  coefficient  (M  unn  1966);  ul  = 


(i.e.,  the  friction  velocity 


squared,  is  a  measure  of  the  surface  stress  associated  with  the  drag);  is  the  mean  wind 
velocity  squared;  and  A  (as  defined  above)  is  the  leaf  area  density  for  the  forest  canopy.  The 
assumption  of  horizontal  homogeneity  also  implies  that  within  the  canopy  itself,  far  enough 
away  from  leading  or  trailing  edge  transition  zones  (i.e.,  in  equilibrium  flow),  there  are  no 
appreciable  changes  in  the  horizontal  pressure  gradient;  in  other  words,  the  mean  horizontal  wind 

dU  d  U  dp 

speeds  do  not  decelerate, - <  0,  or  accelerate, - <0,and  —  is  approximately  constant. 

d  X  d  X  d  X 

Also,  through  scale  analysis,  one  finds  that  the  accelerations  due  to  the  horizontal  pressure 
gradient  forces  are  on  the  order  of  10"''  ms"^  Therefore,  in  the  upper  part  of  the  canopy,  when 
conditions  are  such  that  turbulent  shearing  stresses  and  drag  forces  are  dominant,  the  following 
expression  holds: 


d  u'  w' 
d  z 


\ciA  1/2 


(22) 


Solving  this  complex  equation*  results  in  the  well  known  exponential  or  extinction  profile  for 
wind  speeds  in  the  upper  portion  of  the  forest  canopy,  which  can  be  written  as  described  by 
Inoue  (1963),  Cionco  (1965),  Shinn  (1971),  and  Albini  (1981)  as 


U  =  Uf^  exp  n  1 


z 

~h 


(23) 


By  the  mixing  length  hypothesis,  the  turbulent  shearing  stress  is  related  to  the  mean  velocity  gradient  so  that  the  dependent 

-  2 


variable  on  the  left  hand  side  of  equation  (22)  can  be  rewritten  as  li  w'  =  1‘ 


dU 


dz 


dU  2  T  j  ■ 

- =  /  -  .  Its  derivative  with  respect 

dz  dz 


.  ~,rldU  d^U  1  —2 

to  height  yields  21  — — —C^AU  or 


dz  dz^ 


dU  d'^U 
dz  dz^ 


_  jj  1  Pj-qj^  (Pjg  expression,  one  can  show  that  equation  (23) 


4/ 


is  a  satisfactory  result. 
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1 


where  Uf^  is  the  mean  wind  speed  at  the  forest  top,  n  =  h 


CdA 

4/2 


is  the  extinction  coefficient, 


which  is  a  function  of  leaf  area  density  and  surface  drag,  and  /  is  mixing  length,  often  used  in 


parameter! zations  to  estimate  the  momentum  flux  (i.e.,  m'w'  = 


du 

du 

dz 

dz 

The  extinction 


coefficient,  n  also  has  been  referred  to  as  the  canopy  flow  index  (see  Table  1),  the  values  of 
which  have  been  determined  experimentally  for  various  forest  canopy  types  (Shinn  1971). 
Alternately,  M  eyers  (1987)  and  M  eyers  etal.  (1998)  suggest  the  foil  owing  modification  to 
equation  (23) 


U  =  Uh  exp  1  -  ,  (24) 

h 

P  1 

where  =  LA/ and  =1  ,  where  P  =  1,  2,  or  3  correspond  to  one  of  the  three  general 

profile  shapes  (leaf  area  distributions)  for  forest  canopies,  as  discussed  earlier.  Finally,  the 
extinction  profiles  described  here  could  be  similarly  applied  to  estimate  the  wind-turning  angle 

[i.e.,  o  =  ^  3  ]  from  the  subcanopy  wind  velocity  components  M  =  L'sin  q  and 

u 

v  =  Ucos  0  (Shinn  1971). 


Table  1.  Canopy  flow  Indices  determined 
experimentally  at  various  forest 
sites,  as  reported  by  Shinn  (1971) 


Gum-Maple 

4.42  ±  1.05 

Maple-Fir 

4.03  ±  0.69 

Jungle 

3.84  ±  1.52 

Spruce 

2.74  ±  1.29 

Oak-Gum 

2.68  ±  0.66 

In  the  lower  portion  of  the  forest  canopy  (i.e.,  through  the  open  trunk  spaces  where  air  flow  is 
considerable  less  restricted)  wind  speeds  are  said  to  be  influenced  by  drag  at  the  ground  surface, 
thermal  stability  (inverted  or  convective),  wind  gusts  that  penetrate  through  clearings  at  the 
forest  top;  and  to  a  lesser  extent  (in  a  large  uniform  expanse),  horizontal  advection  (Shinn  1971; 
Shaw  1977;  and  Holland  1989).  In  addition,  a  low-level  wind  maximum  is  often  observed  at  a 
height  around  the  base  of  the  live  branches,  the  magnitude  of  which  is  strongly  influenced  by  the 
forest  type  and  leaf  area  density.  Atthis  level,  the  turbulent  wind-shearing  stress  may  be  depleted 
because  momentum  is  absorbed  most  strongly  through  the  upper  canopy  layers  (i.e.,  within  the 
tree  crown  [Shaw  etal.  1974]).  Therefore,  equation  (21)  may  be  given  instead  as 
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(25) 


0  = 


d  p 

d  X 


-CA 

2  ^ 


where  wind  speeds  will  tend  to  increase  as  the  pressure  drag  forces  decease  across  the  trunk 
space  in  the  absence  of  leaves  and  branches.  It  has  been  suggested  that  the  low-level  wind  speed 
maximum  in  the  lower  portion  of  the  forest  canopy  occurs  in  a  transition  layer  between  the  upper 
canopy— extinction  type  profile  and  the  log-law  type  wind  speed  profile  that  we  would  expect 
extends  from  this  height  to  the  ground. 

3.4  Ist-Order  Closure  and  Higher-Order  Closure  Models  of  the  Forest  Canopy 

Alternately,  several  authors  have  reported  on  the  use  of  first-order  closure  models  (Li  etal.  1990; 
Wilson  etal.  1998;  Wilson  and  Flesch  1999;  and  Pinard  and  Wilson  2000)  and  higher-order 
closure  models,  to  include  large  eddy  simulations,  to  estimate  wind  speeds  and  turbulence  within 
and  above  forest  canopies  (Wilson  and  Shaw  1977;  Y  amada  1982;  M  eyers  and  Paw  U  1986; 

M  eyers  1987;  W  ilson  1988;  Paw  U  and  M  eyers  1989;  Shaw  and  Schumann  1992;  Shen  and 
Leclerc  1997;  KatuI  and  Albertson  1998;  Pyles  etal,  2000;  Albertson  etal.  2001;  and  Li  2001). 
Primarily,  the  governing  equations  for  these  models,  neglecting  coriolis*  and  buoyancy  forces, 
have  been  expressed 


lAi 

dXj 


=  0  , 


(26) 


for  continuity  or  conservation  of  matter,  where  uj  is  the  j-component  of  the  vector  (m,v,vv),  and 
Xj  is  the  j -component  of  the  vector  {x,y,z),  and. 


d  t 


-I- 


—  d  U;  d  UjU  j 

Uj - -  +  - - 

d  Xj  d  Xj 


d  p 
d  Xi 


d  Xi 


dhj 


(27) 


for  the  mean  flow  (e.g.,  Wilson  and  Shaw  1977),  where  the  overbar  and  primed  variables 
indicate  the  mean  and  fluctuating  components  of  the  given  quantity.  Equation  (27)  contains  the 
following  terms  from  I  eft  to  right:  1)  time  tendency  of  the  mean  wind,  2)  advection  of  the  mean 
wind,  3)  eddy  flux  divergence,  4)  the  mean  pressure  gradient,  5)  changes  in  the  fluctuating 
component  of  static  pressure  (i.e.,  the  pressure  drag  force),  6)  the  viscous  drag  accelerations  on 
the  mean  flow,  and  7)  viscous  drag  accelerations  on  the  canopy.  The  5th  and  7th  terms  represent 
"the  summation  of  all  the  differential  pressure  forces  (and  velocity  gradients)  on  the  upstream 
and  downstream  sides  of  the  elements  of  the  (forest)  canopy  (Wilson  and  Shaw  1977)."  Where 
pressure  forces  are  considered  the  main  portion  of  the  total  drag  from  the  forest  canopy,  then  the 
following  closure  assumption  is  made 


In  meteorology,  the  coriolis  force  is  an  apparent  (deflecting)  force  on  moving  air  parcels  in  the  atmosphere  that  arises  due  to 
the  earth’s  rotation  so  that/=  -2  V,  which  is  to  the  right  of  motion  in  the  northern  hemisphere  (Huschke  1959).  Here  is  the 
angular  velocity  of  the  earth  and  V  is  velocity.  Other  than  a  few  authors  (e.g.,  Shinn  1971,  Holland  1989,  and  Wilson  and  Flesch 
1999),  most  exclude  the  effect  of  the  coriolis  force  as  having  negligible  effect  on  the  scales  of  motion  considered  (i.e.,  windflow 
through  the  forest  subcanopy  layer  see  equation  [25]). 
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(28) 


In  addition,  the  most  common  first-order  turbulence  closure  for  the  momentum  flux  divergence 
(term  3)  in  equation  (27)  can  be  expressed,  as  described  by  Li  etal.  (1990)  as 


d  UiUj  3  ^  ^Uj  ^ 

3  Xj  dxj  ^  dxj  dxj 


(29) 


where  Kaj)  is  the  viscosity  (or  eddy  exchange  coefficient)  given  by 

7  I  /' 

K(..^  =  r  — -  +  — ^  ,  and  /  is  a  characteristic  mixing  length,  which  has  been 

y  dxj  dxi 

parameterized  in  various  ways  as  a  function  of  the  forest  height,  canopy  density,  zero-plane 
displacement,  and  roughness  (Li  etal.  1990;  Wilson  etal.  1998;  and  Pinard  and  Wilson,  2001). 
Alternately,  Wilson  etal.  (1998)  specify  the  eddy  exchange  coefficient  in  terms  of  a 

characteristic  mixing  length  scale,  \x  I  ^  and  the  square  root  of  the  turbulent  kinetic 

energy,  where  =  ^  .  As  a  result, 

2  i 

K=  mix  ,  (30) 


2 

where - is  the  ratio  of  the  equilibrium  shearing  stress  and  turbulent  kinetic  energy  at  the 

tke{h) 

forest  canopy  top. 


4.  Vertical  Profiles  of  Temperature  in  Forests 


The  forest  canopy  is  an  active  surface  for  heat  and  moisture  exchanges,  and  therefore,  like  over 
open  fields,  temperature  profiles  above  the  canopy  top  will  form  with  inversions  at  night  and 
unstable  lapse  rate  conditions  during  the  daytime  (M  unn  1966).  However,  temperature  profiles 
above  forests  will  not  normally  maintain  large  (positive  or  negative)  gradients  because  of 
downward  mixing  (see  discussion  on  wind  speeds  above  the  canopy  top  in  section  3.2).  Also, 
inversion  conditions  will  sometimes  develop  during  the  daytime  above  forests  due  to  evaporative 
cooling  from  leaf  transpiration.  In  contrast,  within  the  forest  canopy,  one  mostly  finds  inversion 
conditions  during  the  day  and  isothermal  or  slightly  unstable  lapse  conditions  at  night  (Figure  5). 
Bergen  (1974)  describes  the  diurnal  variation  of  temperatures  within  the  forest  canopy  in  the 
following  manner:  a  local  temperature  maximum  develops  between  0900  and  1000  LT,  which 
moves  downward  into  the  canopy  as  solar  altitude  increases;  the  maximum  penetration  of  heat  is 
reached  near  1200  LT  at  a  height  of  about  0.7  A;  variation  in  temperatures  through  the  sub- 
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Figure  5.  Observed  mean  temperatures  (C)  through  a  tropical  rain  forest 
canopy  for  different  times  of  day  (based  on  Bayton  1963). 


canopy  space  is  complex,  but  inversion  conditions  are  often  observed  throughout  the  daytime 
period. 

Likewise  Raynor  (1971)  and  Hoskeretal.  (1974)  describe  the  forest  canopy  for  heat  absorption 
and  radiative  exchanges.  During  the  daytime,  the  canopy  captures  incoming  solar  radiation  and  a 
layer  of  warm  air  develops  that  propagates  downward  through  normal  mixing  processes.  As  an 
example.  Figure  6  shows  a  typical,  midday  temperature  profile  for  a  pine  forest  stand  on  a  clear 
day.  M  aximum  temperatures  are  shown  to  coincide  with  the  height  of  maximum  leaf  density 
(and  minimum  wind  speeds;  see  Figure  2),  whereas  a  local  minimum  in  temperature  occurs 
through  the  trunk  space  below  the  live  branches.  Under  cloudy  conditions  with  and  without 
rainfall,  Floskeretal.  (1974)  observed  mostly  isothermal  conditions  within  a  pine  forest  canopy 
throughout  the  day. 
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Figure  6.  Typical  vertical  temperature  profile  structure  within 

and  above  a  pine  forest  stand  (based  on  Bergen  1974). 


In  contrast  at  night,  the  forest  canopy  loses  heat  and  is  cooled  further  (from  the  upper  crown  area 
downward)  by  subsidence  so  that  isothermal  conditions  or  weak  negative  lapse  rate  conditions 
prevail.  This  was  shown  in  the  time  sequence  of  vertical  profiles  reported  by  Huisman  and 
Attenborough  (1991).  A  second  temperature  inversion  is  likely  to  form  close  to  the  ground  at 
night  due  to  surface  radiative  cooling  of  the  soil,  as  had  been  observed  by  Fontan  etal.  (1992) 
and  Lee  etal.  (1997). 

The  conservation  of  energy  (or  enthalpy,  Cp~)  equation  can  be  written  as  described  by  Businger 
(1982)  and  Garratt(1992), 

d~  -d~  ^“7  1  dRj 

- h  U  j - — - h  Y - T  -  (31) 

dt  ■'  dxj  dxj  dxj  Cp  dxj 

where  the  terms  from  I  eft  to  right  are  1)  the  time  tendency  of  the  mean  potential  temperature  ( ~), 
2)  advection  of  potential  temperature  by  the  mean  wind,  3)  eddy  heat  flux  divergence,  4)  the 
viscous  dissipation  of  heat  and  5)  the  net  radiative  flux  divergence  through  the  forest  canopy. 

The  heat  flux  divergence  term  can  be  expanded  as 


d  u, 


d  u' 


d  V 


d  w' 


dx; 


dx 


■  +  ■ 


dy 


■  +  ■ 


dz 


(32) 
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Bergen  (1974)  suggested  that 


d  w' 

[I  n  (t^  Tf)  A  ,  (33) 

where  n  is  the  heat  transfer  coefficient  for  the  leaf  (pine  needle)  surface,  Ta  is  the  air  temperature, 
T/ is  the  leaf  (needle)  surface  temperature,  and  A  is  the  foil  age  surface  area  in  the  volume. 

M  oreover,  Bergen  (1974)  indicated  that  the  value  for  (r^-T/)  would  depend  on  the  solar  altitude, 
the  ratio  of  diffuse  to  direct  solar  radiation,  and  cloud  cover.  In  addition,  the  angle  between  the 
wind  velocity  at  the  canopy  top  and  the  solar  azimuth  would  be  an  important  factor  because  of 
the  influence  (on  temperature)  from  ventilation  through  opening  at  the  canopy  top  and  additional 
heating  in  direct  sunlight.  The  reverse  situation  is  considered  for  decreased  wind  flow  and 
shaded  areas  in  the  lee  of  dense  crowns.  Intermittent  wind  gusts  through  openings  at  the  top  of 
the  forest  canopy  are  also  important  because  they  will  bring  downdrafts  of  cooler  and  drier  air 
from  above  into  the  canopy  (Denmead  and  Bradley  1985). 


5.  Acoustics  in  Forests 


5.1  Approximation  of  Sound  Speed  Profiles  in  Forests 

Figure  7  shows  derived  sound  speed  profiles  within  and  above  a  10  m  (pine  forest)  canopy  given 
the  characteristic  variations  in  wind  speed  and  temperature  that  were  described  earlier  in  sections 

3.1  and  4.  Using  equations  (10-12)  from  section  1,  the  three  plotted  profiles  represent  from  left 
to  right  1)  the  effective  sound  speed,  along  the  path  of  propagation  from  the  source  to  the 
receiver  (i.e.,  upwind  propagation  [see  Figure  1]),  where  r  =  90°  is  the  bearing  of  the  receiver 

from  the  source,  and  w  =  90°  is  the  direction  of  the  wind  from  North;  2)  the  speed  of  sound  in 
air,  CO,  from  equation  (10),  which  is  the  speed  of  sound  in  the  absence  of  wind  or  crosswind 
propagation,  and;  3)  the  effective  sound  speed,  Cejf,  where  r  =  90°  and  ^  =  270°  (i.e., 

downwind  propagation).  The  effective  speed  of  sound  for  downwind  propagation,  =  270°,  is 
of  a  greater  magnitude  across  the  entire  profile  than  in  the  opposite  case  (i.e.,  where  r  =  90°  and 

qw  =90°).  Flowever,  the  sound  speed  profile  here  is  flatter  (and  its  slope  above  the  canopy  top  is 
positive)  in  comparison  to  the  other  profiles.  The  influence  these  profile  variations  in  sound 
speed  have  on  short-range  acoustic  attenuation  will  be  discussed  briefly  in  the  following  section. 
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Figure  7.  Derived  speed  of  sound  profiles  within  and  above  a  10-m  pine 
forest  canopy  given  the  characteristic  variations  in  wind  speed 
and  temperature  that  were  shown  in  Figures  (2)  and  (6).  The  three 
plotted  profiles  are  described  in  the  text. 


5.2  Approximation  of  Short-Range  Acoustic  Attenuation  in  Forests 

Attenuation  of  sound  waves  in  forests  is  said  to  involve  three  main  phenomena:  1)  interference 
between  direct  and  ground  reflected  acoustic  waves;  2)  scattering  by  tree  trunks  and  branches, 
the  ground,  and  turbulence;  and  3)  absorption  by  the  trees,  leaves,  branches,  the  ground,  and  the 
air  (Erying  1946;  Embleton  1963;  Aylor  1972;  Burns  1979;  Fricke  1984;  Attenborough  1985; 
Price  etal.  1988;  Huismann  and  Attenborough  1991;  and  Sakai  etal.  1998,  2001).  Ground 
impedance,  which  affects  sound-wave  reflection  and  absorption,  has  been  found  to  be  a  function 
of  flow  resistivity  and  ground  surface  porosity.  Ground  impedance  has  been  reported  to  affect 
sound  attenuation  mainly  in  the  250-500  Hz  frequency  range  (Fricke  1984;  Price  etal.  1988;  and 
Sakai  etal.  2001).  Flow  resistivities  for  ground  surfaces  under  typical  forest  stands  have  been 
reported  as  =  20-80  kPa  s  m"^  which  are  much  lower  in  comparison  to  those  reported  for  grass 

(  =  150-300  kPa  s  m"^)  or  soil  covered  surfaces  (  =  300-800  kPa  s  m"^)  (Embleton  etal., 
1983).  Likewise,  the  forest  floor,  which  consists  of  humus,  leaves  (needles),  and  other  biological 
litter,  has  been  reported  to  have  values  for  porosity  of  about  =  0.825,  which  are  considerably 

higher  than  those  reported  for  grass  fields  (  =  0.675)  or  unpacked  soils  (  =  0.575) 
(Attenborough  1985). 

At  higher  frequencies  (/>  1  to  2  kHz),  sound  attenuation  in  forests  has  been  observed  as  being 
mainly  due  to  scattering  by  tree  trunks  and  branches  and  absorption  by  foliage  (Embleton  1963; 
Aylor  1972;  Burns  1979;  Fricke  1984;  Price  etal.  1988;  and  Huismann  and  Attenborough  1991). 
These  effects  on  sound  attenuation  at  higher  frequencies  may  be  as  important  as  refraction 
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effects  due  to  temperature  inversions  and  wind  speeds  within  the  forest  canopy  (Huismann  and 
Attenborough  1991). 

To  provide  an  example,  we  have  used  the  characteristic  profiles  of  wind  speed  and  temperature, 
which  were  shown  in  Figures  (2)  and  (6)  as  input  to  a  flat-earth,  nonturbulent,  acoustic  propagation 
model  called  WSCAFFIP  (the  Windows  (version)  Scanning  Fast  Field  Program).  WSCAFFIP  is  a 
numerical  code  developed  for  assessing  environmental  effects  on  short-range  acoustic  attenuation 
(Noble  and  Marlin  1995).  WSCAFFIP  determines  acoustic  attenuation  as  relative  sound  pressure 
loss  with  range  and  azimuth  for  a  given  frequency  and  source- receiver  geometry.  WSCA  FFI P 
contains  propagation  algorithms  to  represent  the  effects  of  atmospheric  refraction,  diffraction, 
absorption,  and  reflection  (ground  impedance)  on  acoustic  transmission.  Table  2  lists  the  model 
parameters  for  an  initial  approximation  of  short-range  acoustic  attenuation  in  forests.  The  distance 
of  the  receiver  from  the  source,  source  height,  and  receiver  height  were  500  m,  2  m,  and  1  m, 
respectively.  We  have  calculated  relative  attenuation  at  50, 100,  500,  and  1000  Hz  for  upwind 
propagation  (  «  =  90°and  =90°)  and  downwind  propagation  (  «  =  90°and  w  =  270°).  Figures 
shows  the  WSCAFFIP  results  for  the  sound  speed  profiles  described  in  section  5.1. 


Table  2.  WSCAFFIP  model  parameters  for  approximation  of  short-range  acoustic 
attenuation  In  forests. 


Parameter 

Symbol 

Value 

Distance  of  receiver  from  source 

Jt 

500  m 

Range  resolution 

;v: 

10  m 

Bearing  of  receiver  from  source 

R 

90° 

Source  height  above  ground 

hs 

2  m 

Receiver  height  above  ground 

hn 

1  m 

Frequency  of  interest 

f 

50  Hz;  100  Hz,  500  Hz;  1000  Hz 

Ground  porosity 

0.850 

Flow  resisitivity 

23  kPa  s  mr^ 

Bearing  of  the  wind  from  North 

W 

90° ;  270° 
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Figure  8.  WSCAFFIP  numerical  approximations  of  short-range  acoustic 

attenuation  within  a  continuous  forest  stand.  Calculations  of  relative 
attenuation  are  shown  at  50, 100,  500,  and  1000  FI z  for  r  =  902and 
„  =  902  and  2702. 


Acoustic  waves  within  and  above  the  forest  canopy  will  tend  to  be  refracted  upward  as  the 
effective  sound  speed  decreases  with  height  and  refracted  downward  as  sound  speed  increases 
with  height.  Therefore,  we  might  expect  the  amounts  of  attenuation  to  depend  on  the  strength  and 
locations  of  sound  speed  profile  inversions  between  the  ground  and  the  forest  canopy  top.  Because 
the  sound  speed  profile  for  =  270s  was  comparatively  flatter  than  that  for  „  =  90s  we  would 

expect  to  see  slightly  greater  attenuation  with  range  for  =  270s  at  each  frequency  as  a  result. 

Figure  8  shows  this  result  satisfactorily,  even  though  ground  impedance,  more  that  sound  speed 
profile  shape,  may  have  affected  the  model  results  at  frequencies  <  500  Hz  (Embleton  1996). 

Finally,  in  preparing  this  reportfor  publication  a  reviewer  offered  the  following  comments:  The 
behavior  of  the  results  shown  in  Figure  8  can  be  roughly  understood  in  terms  of  'ducting'  of 
acoustic  modes.  An  acoustic  duct  exists  between  the  surface  and  about  7  m.  At  lower 
frequencies,  50  and  100  Hz,  acoustic  waves  may  be  too  long  for  any  significant  ducting  to  occur. 
AtSOO  Hz,  there  appears  to  be  a  single  trapped  mode;  whileatlOOO  Hz,  itappears  that  two 
modes  are  trapped,  creating  the  i  nterference  pattern  (Wilson  2002).  Thus,  it  is  the  'ducting'  of 
acoustic  modes  inside  and  above  the  forest  canopy  that  is  greatly  influenced  by  variations  in 
local  wind  speed  and  temperature  profile  structure. 


19 


6.  Summary  and  Conclusions 


We  have  initiated  research  to  examine  the  calculation  of  local  speed  of  sound  in  the  forest 
environment  for  future  military  applications  of  acoustic  information  inforests.  Wegathered 
information  from  the  literature  on  measured  and  modeled  micrometeorology  in  forests  to  include 
reports  on  wind  speed  and  temperature  profile  structure,  turbulence,  leaf  area  distributions,  and 
estimates  of  surface  roughness  and  displacement.  (An  Appendix  summarizes  the  results  of  the 
literature  survey  from  which  the  report  was  produced.)  We  also  derived  initial  approximations 
for  sound  speed  and  short-range  acoustic  attenuation  through  a  continuous  forest  stand  based  on 
two  characteristic  profiles  for  wind  speed  and  temperature  inside  and  above  the  forest  canopy. 

As  discussed  in  section  5.2,  attenuation  of  sound  waves  in  forests  is  affected  by  interference 
between  direct  and  ground  reflected  acoustic  waves,  refraction,  scattering,  and  absorption. 
Because  acoustic  waves  tend  to  be  refracted  upward  as  the  effective  sound  speed  decreases  with 
height  and  refracted  downward  as  sound  speed  increases  with  height,  we  expect  that  attenuation 
will  depend  greatly  on  the  strength  and  locations  of  sound  speed  profile  inversions  between  the 
ground  and  the  forest  canopy  top.  Several  initial  model  calculations  showed  this  result  satisfactorily. 
Weexpectthattheinformationgatheredasaresultof  thissurvey  will  be  useful  to  helpevaluate 
various  meteorological  computer  models  and  techniques  for  determining  vertical  profiles  of  sound 
speed  through  the  atmosphere  in  forests. 
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Appendix  -  Literature  Survey 


This  Appendix  contains  several  tables  to  list  the  numerous  citations  associated  with  the  topics 
discussed  in  the  report,  i.e., 

•  Battlefield  acoustic  models  and  acoustic  sensor  arrays, 

•  Speed  of  sound  physics, 

•  Observed  micrometeorology  in  forests, 

•  Computer  models  to  provide  meteorological  profile  and  turbulence  information  within  and 
above  forests, 

•  Forest  characteristics  to  include  leaf  area  density,  roughness,  and  displacement  height,  and 

•  Measurements  and  approximations  of  short-range  acoustic  attenuation  in  forests. 
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Table  A- 1.  The  Future  Combat  System's  battlefield. 


U.S.  Army  Acoustic  Systems; 

a.  Surveillance, 

b.  Tracking, 

c.  Detection, 

d.  Identification, 

e.  Classification, 

f.  Remote  sentry,  and 

g.  Intelligence  gathering  of  military  targets. 

New  U.S.  Army  acoustic  systems  include  unattended 
ground  sensors  (UGS)  in  remote  netted  acoustic/seismic 
sensor  arrays  (RNADS). 

Benn  and  Hilt  (1975); 

Hafner  et  al.  (1976); 

Zappi  (1978);  Pilette  et  al.  (1979); 
Thompson  and  Scherer  (1991); 

Fong  and  Srour  (1994); 

Rosenthal  and  Stevens  (1994); 

Carnes  and  Morgan  (1995); 

Srour  and  Robertson  (1995); 

Wilson  (1997);  Eom  et  al.  (1999); 

Fopez  et  al.  (1999);  Wellman  (1999); 
Wellmann  and  Srour  (1999); 

Young  et  al.  (1999); 

Becker  and  Gudesen  (2000); 

Depireux  and  Shamma  (2000); 

Mays  and  Price  (2000); 

Mays  and  Vu  (2000). 

Army  Acoustic  Computer  Models; 

Smith  (1989); 

Sparrow  (1991); 

Goal;  To  determine  point-to-point  losses  in  acoustic 

West  et  al.  (1991); 

transmission; 

West  et  al.  (1992); 

Auvermann  and  Goedecke  (1993); 

a.  Spherical  spreading  model 

Wilson  (1993);  Noble  (1994); 

b.  Parabolic  Equation  (PE)  model 

Auvermann  et  al.  (1995); 

c.  Fast  field  Program  (PEP) 

Noble  and  Marlin  (1995); 

d.  Scanning  Fast  Field  Program  (SCAFFIP) 

Kotiaho  (1996); 

e  .  Acoustic  Multistream  Propagation  Program 

Goedecke  et  al.  (1997); 

(AMPP) 

Wilson  (1998ab); 

f.  Acoustic  Battlefield  Decision  Aid  (ABFA) 

Goedecke  et  al.  (2000); 

g.  Turbulence  spectral  model 

Wilson  (2000); 

Wilson  and  Szeto  (2000); 

Atmospheric  turbulence  and  refraction  effects  influence  the 

Auvermann  (2001). 

retrieval  and  interpretation  of  acoustic  signals. 
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Table  A -2.  The  atmospheric  physics  of  sound  speed  (based  on  Fleagle  and  Businger  1963). 


Change  in  pressure  {dp)  over  a  unit  cross-sectional  area  of  a 
fluid  containing  a  sound  wave  (i.e.,  between  the  plane  of 
compression  and  plane  of  rarefaction); 

dc 

(A-1)  dp  =  ds  — 

dt 

Ratio  of  pressure  change  to  volume  change  for  an 
adiabatic  (or  isentropic)  process; 

(A-2)^= 

d  c„ 

Equation  of  state  (EOS)  for  an  ideal  gas; 

(A-3)  =  — , 

M 

where  R  =  8314.32  J  is  the  universal  gas  constant  and 

M  is  molecular  mass. 

Speed  of  sound  in  air; 

.X  1  RT 

(A-4)  c  =  , 

V  M 

(A-5)  =  — —  1.403  is  the  ratio  of  specific  heats, 
r'v 

(A-6)  c„  =  331.29  ms  '  is  the  reference  value  for 
speed  of  sound  for  standard  dry  air  at  0  °C  and  a 
barometric  pressure  of  1013.25  mb  (Wong  1986). 

Ratio  of  specific  heats  and  molar  mass; 

Effective  sound  speed  (Figure  1); 

(A-7) —  =  0.04833  + (a  0.023)4,’ 

M 

(A-8)  =  C„  +  M  cos  (  ^  , 

where  accounts  for  increases  (decreases)  of  sound 

A,  =9.2  X  10  ^  +  5.5  jc  10  +  4.25  x  10 

speed  due  to  variations  in  wind  velocity  (Noble  and 

h  is  humidity  (in  the  range  0.0  to  1 .0),  and  f  is  air  temperature 

Marlin  1995  and  Osteshev  1997) 

in  degrees  Celsius  (Wong  and  Embleton  1984,  1985) 
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Table  A -3.  Vertical  profiles  of  wind  speed  in  forests. 


The  study  of  winds  in  forests: 

A.  To  better  understand  the  surface  layer  exchanges  of  heat,  moisture  (water  vapor),  and  carbon  dioxide,  which 
not  only  influence  regional  and  local  climate,  but  also  affect  the  forest  ecosystem  for  growth  and  regeneration; 

B.  For  the  study  of  transport,  diffusion,  and  deposition  of  air-borne  pollutants,  trace  gases,  and  aerial  sprays; 

C.  To  assess  the  impact  of  wind  damage  on  trees  and  tree  harvests  by  severe  storms. 


Survey  of  the  literature: 

A.  Meteorological  measurements; 

B.  Model  forests  in  wind  tunnels; 

C.  Atmospheric  computer  models  of  varying  complexity. 


Allen  (1968);  Bergen  (1971);  Martin  (1971); 

Oliver  (1971);  Raynor  (1971);  Denmead  and  Bradley  (1985); 
Amiro  and  Davis  (1988);  Baldocchi  and  Hutchinson  (1988); 
Baldocchi  and  Meyers  (1988);  Lee  and  Black  (1995). 

Meroney  (1968);  Sadeh  et  al.  (1971);  Chen  et  al.  (1995). 


Monin  and  Obukhov  (1954);  DeMarrais  (1959); 
Inoue  (1963);  Cionco  (1965);  Barr  (1971);  Shinn  (1971); 
Busch  (1973);  Dyer  (1974);  Benoit  (1977);  Shaw  (1977); 
Touma  (1977);  Wilson  and  Shaw  (1977);  Irwin  (1979); 
Garratt  (1980,  1992);  Albini  (1981);  Shaw  and 

Peri  era  (1982);  Yamada  (1982);  Meyers  and  Paw  U  (1986); 
Meyers  (1987);  Wilson  (1988);  Paw  U  and  Meyers  (1989); 
Diet  al.  (1990);  Shaw  and  Schumann  (1992);  Kaimal  and 
Finnegan  (1994);  Zoumakis  (1993,  1994);  Finnegan  and 
Brunet  (1995);  Lo  (1995);  Massman  (1982);  Shen  and 
Leclerc  (1997);  Katul  and  Albertson  (1998);  Meyers  et  al. 
(1998);  Wilson  et  al.  (1998);  Wilson  and  Flesch  (1999); 
Pyles  et  al.  (2000);  Albertson  et  al.  (2001);  Li  (2001); 
Pinard  and  Wilson  (2001). 


Main  features  of  wind  speed  profile  structure  (Figure  2): 

1)  A  low-level  wind  maximum  at  -0.1  to  03h,  where  h  is  the  height  of  the  canopy  top; 

2)  A  layer  of  minimum  wind  speed  at  -0.6  to  0.8/i,  which  coincides  with  the  region  of  maximum 
leaf  area  density; 

3)  A  diurnal  time  maximum  of  wind  speed,  which  extends  through  the  entire  canopy  in  the  early  afternoon. 
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Table  A -4.  Wind  speeds  above  a  uniform  forest  stand. 


Equations  for  (idealized)  steady  state, 
homogeneous,  equilibrium  conditions 


Annotation 


Reference 


A.  Wind  flow  with  constant  shearing  stress: 
d  U  M* 


(A-9) 


k(z  d) 


(A-10)  [/  =  —  In  — 

k  zq 


B.  Power-law  exponent  model 

(A-11)  [7,  =c7, 

zi  d 


U  is  the  mean  wind  speed; 

Ut  =  is  the  friction  velocity; 

—  =  u'w'  is  the  shearing  stress; 
is  air  density; 

k  is  Karman’s  constant,  k  0.4  ; 
d  is  zero-plane  displacement; 

Zq  is  surface  roughness  length; 

^  is  a  function  of  stability;  and 
„  is  the  diabatic  influence  function. 


U  and  U  are  the  wind  speeds  at 
heights  Z[  and  Z2  above  h  ,  the 
canopy  top  and  the  exponent  p  is  a 
function  of  wind  speed  and  stability. 


Monin  and  Obukhov  (1954); 
Panofsky  (1963); 

Paulson  (1970); 

Webb  (1970); 

Busch  (1973); 

Dyer  (1974); 

Benoit  (1977). 


DeMarrais  (1959); 

Martin  (1971)  Touma  (1977); 
Irwin  (1979). 


C.  Windflow  within  the  transition  layer  (i.e., 
H  z<  z,): 


(A-12) 


a  u 


3  z  k{z  d) 


Within  the  transition  layer,  additional 
momentum  is  transferred  to  the  winds  due  to 
surface  (canopy-top)  generated  eddies, 
which  effectively  decreased  the  vertical 
profile  gradient  (Figure  4). 

z.  is  the  depth  of  the  transition  layer. 


Garratt(1980,  1992); 
Zoumakis  (1993,1994); 
Kaimal  and  Finnegan  (1994); 
Arya  (2001). 


(A-13)  „=exp  0.7  1  — 

z* 


.  .  ,  1  ^  d  * 

(A-14)  [/  =  —  In  -  m  +  m 

k  Zo 


z. 


3  ,  where  i 
z .  =  1.5/i  to  2.5h, 
roughness  density. 


:  tree  spacing,  or 
depending  on  the 


„  is  defined  in  the  same  manner  as 
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Table  A -5.  Wind  speed  profile  structure  within  a  uniform  forest  stand. 


Equations  for  (idealized)  steady  state, 
homogeneous,  equilibrium  conditions 


Annotation 


Reference 


From  conservation  of  momentum: 


(A-15)  ^  U\ 

d  z  d  X  2 


a.  Momentum  flux  divergence, 

b.  Horizontal  pressure  gradient,  and 

c.  Drag  due  to  friction  imparted  by  the 
leaves  and  branches. 


This  yields. 


(A-16)  ^  =  -QA 

d  z  ^ 


2^  2, 

C  j  =  is  the  surface  drag  coefficient 


A  is  the  leaf  area  density  per  unit  surface  area  of 
ground  below. 


LAI  =  A(z)  dz  is  leaf  area  index;  Generally, 
0 

LAI  =  1  to  5  for  forests. 


Under  conditions  that  the  accelerations  due  to 
turbulent  shearing  stresses  and  drag  forces  in  the 
upper  part  of  the  canopy  are  dominant  (e.g.,  far 
enough  away  from  leading  or  trailing  edge 
transition  zones)  then  equation  (A-16)  holds. 


Barr  (1971); 

Shaw  (1977); 

Wilson  and  Shaw  (1977). 


Pielke  (1984); 

Kaimal  and  Finnigan 
(1994);  Finnigan  and 
Brunet  (1995). 


Upon  integration,  equation  (A-16)  yields, 


(A-17)  (7  =  [/;,  exp  Ml  - 

h 


U the  mean  wind  speed  at  the  forest  top. 


n  =  h 


CjA 

4P 


is  the  extinction  coefficient  or 


Inoue  (1963); 
Cionco  (1965); 
Shinn  (1971); 
Albini  (1981). 


canopy  flow  index,  and  /  is  mixing  length,  i.e.. 


u' w' =  1‘ 


dU 


dz 


dU 

dz 


Alternately, 


(A-18)  U  =  Uh  exp  1  - 

h 


=  LAI, 


=  1 


P  1 


,  and  R=  1,  2,  or  3 


Paw  U  and  Meyers 
(1989); 

Meyers  et  al.  (1998). 


correspond  to  one  of  the  three  general  profile 
shapes  (leaf  area  distributions)  for  forest  canopies 

(Figure  3). 


The  subcanopy  wind  velocity  components: 
(A-19)  u  =U  sin  q 
(A-20)  V  =  U  cos  0 


1  V 

Q  =  tan  —  is  the  subcanopy 
u 

wind-turning  angle 


Shinn  (1971). 
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Table  A -6.  Zero-Plane  displacement,  roughness  length,  and  subcanopy  wind  maxima. 


Equations  for  (idealized)  steady  state, 
homogeneous,  equilibrium  conditions 

Annotation 

Reference 

Minimum  wind  speeds  through  the  tree  crown 
have  been  associated  with  the  zero-plane 
displacement  height. 

LAI  is  the  leaf  area  index  and  P  =  1,  2,  or  3 
correspond  to  one  of  the  three  general  profile 
shapes  for  forest  canopies. 

Shaw  and  Periera  (1982); 
Meyers  (1987); 

Paw  U  and  Meyers  (1989); 
Meyers  et  al.  (1998). 

(A-21)  J  =  h  0.05 -f 

2  20 

From  various  experiments  in  pine  forests, 
d  ~  0.70  -  0.76/1. 

Bergen  (1971); 

Raynor  (1971); 

Oliver  (1971); 

Lee  and  Black  (1993). 

Surface  (aerodynamic)  roughness  (in  equation  A- 
10)  has  been  parameterized  in  the  following 
manner: 

/4/0'25  (p  jt 

(A-22)  =  h  0.23 

lU  0/ 

Roughness  (z„)  varies  in  the  following  manner: 

~  10"^  m  over  snow,  sand,  dry  lakebeds; 

~  10"^  m  over  soils; 

~  0. 1  m  over  farmlands,  tall  grass,  or  shrubs; 

~  1.0  m  for  orchards  and  forests. 

Oke(1978); 

Pielke(1984); 

Rachele  and  Tunick  (1994). 

Through  the  open  trunk  spaces  where  air  flow  is  considerable  less  restricted,  equation  (A- 15)  yields 

(A-23)  0  =  -CjA 

d  X  2  “ 

Shinn  (1971); 

Shaw  (1977); 

Shaw  et  al.  (1974); 

Holland  (1989). 

1 )  Wind  speeds  are  mainly  influenced  by  drag  at  the  ground  surface,  thermal  stability,  wind 
gusts  that  penetrate  through  clearings  at  the  forest  top,  and  to  a  lesser  extent  (in  a  large 
uniform  expanse),  horizontal  advection; 

2)  A  low-level  wind  maximum  is  often  observed  at  a  height  around  the  base  of  the  live  branches, 
the  magnitude  of  which  is  strongly  influenced  by  the  forest  type  and  leaf  area  density; 

3)  The  turbulent  wind  shearing  stress  may  be  depleted  because  momentum  is  absorbed  most 
strongly  through  the  upper  canopy  layers  (i.e.,  within  the  crown); 

4)  As  a  result,  the  low-level  wind  speed  maximum  occurs  in  a  transition  layer  between  the 
upper-canopy  extinction  type  profile  and  the  log-law  type  profile  that  we  would  expect 
extends  from  this  height  to  the  ground. 
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Table  A -7.  Ist-Order  closure  models  and  higher-order  closure  models  to  Include  large  eddy  simulations. 


Equations 


Continuity  (Conservation  of  Matter); 


(A-24) 


d  Uj 


0 


Equation  for  the  Mean  Elow; 


(A-25) 


a  t 


-  a  U: 

+  U:  - - 


a  UjUj 


dp  dp 

a  X,.  a  X,. 


a.  Time  tendency  of  the  winds, 

b.  Advection  of  the  mean  wind, 

c.  Momentum  flux  divergence, 

d.  Mean  pressure  gradient, 

e.  Pressure  drag  force, 

f  Viscous  drag  force  on  the  mean  flow,  and 
g.  Viscous  drag  force  from  the  canopy  elements. 


Reference 


a^M;  a^M; 


Wilson  and  Shaw  (1977); 
Yamada(1982); 

Wilson  (1988); 

Meyers  and  Paw  U  (1986); 
Meyers  (1987); 

Paw  U  and  Meyers  (1989); 
Shaw  and  Schumann  (1992); 
Shen  and  Leclerc  (1997); 
Katul  and  Albertson  (1998); 
Pyles  et  al.  (2000); 

Albertson  et  al.  (2001); 

Li  (2001). 


When  the  mean  pressure  gradient  is  considered  negligible,  and  pressure  drag 
forces  are  considered  the  main  portion  of  the  total  drag  from  the  forest  canopy, 
the  following  closure  assumption  is  often  used 


(A-26) 


d  P 

a  X; 


-  CjA  u} 


First-order  turbulence  closures  for  the  momentum  flux  divergence; 


(A-27) 


a  UM  : 

'  J 

a  xj 


where 

(A-28) 

Alternately, 

(A-29)  K  =  ,,,, 


dxj  dxj  ax,. 


dxj  ax,. 


2 

u* 

tke(h) 


tke{z) 


1 

7 


Li  etal.  (1990); 

Wilson  et  al.  (1998); 
Wilson  and  Flesch  (1999); 
Pinard  and  Wilson  (2001). 


where  tke  =  ^  u]  and 


u* 


tke(h) 


is  the  ratio  of  the  equilibrium  shearing  stress 


and  turbulent  kinetic  energy  at  the  forest  canopy  top. 
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Table  A -8.  Vertical  profiles  of  temperature  in  forests. 


Equations 

Reference 

The  conservation  of  energy  (or  enthalpy,  )  equation; 

Businger  (1982); 

Garratt  (1992). 

onx  3“  —a”  ^  1  ^Rj 

(A-30)  u  j  ^  T  ^ 

dt  ■'  dxj  dxj  dxj  Cp  dxj 

a.  Time  tendency  of  the  mean  potential  temperature  (  ); 

b.  Advection  of  potential  temperature  by  the  mean  wind; 

c.  Eddy  heat  flux  divergence; 

d.  Viscous  dissipation  of  heat; 

e.  Net  radiative  flux  divergence  through  the  forest  canopy. 

a  w' 

(A-31)  — -  It  n  Tf)A 

dz 

where 

n  is  the  heat  transfer  coefficient  for  the  leaf  (pine  needle)  surface, 
is  the  air  temperature, 

Tf  is  the  leaf  (needle)  surface  temperature,  and 

A  is  the  foliage  surface  area  in  the  volume. 

Bayton  (1963); 

Munn  (1966); 

Raynor  (1971); 

Bergen  (1974); 

Hosker  et  al.  (1974) 

Denmead  and  Bradley  (1985); 

Huisman  and  Attenborough  (1991); 
Eontan  et  al.  (1992); 

Lee  etal.  (1997). 

Tj)  will  depend  on; 

a.  Solar  altitude; 

b.  Ratio  of  diffuse  to  direct  solar  radiation; 

c.  Cloud  cover;  and 

d.  Angle  between  the  wind  velocity  at  the  canopy  top  and  the  solar 
azimuth. 

Main  features  of  temperature  profile  structure  within  and  above  forests,  Figures  (5)  and  (6); 

1)  Inversion  conditions  during  the  day  and  isothermal  or  slightly  unstable  lapse  conditions  at  night; 

2)  Inversion  conditions  will  sometimes  develop  during  the  daytime  above  forests  due  to  evaporative  cooling 
from  leaf  transpiration; 

3)  Temperature  profiles  above  forests  will  not  normally  maintain  large  (positive  or  negative)  gradients; 

4)  The  forest  canopy  captures  incoming  solar  radiation  and  a  layer  of  warm  air  develops  that  propagates 
downward  through  normal  mixing  processes; 

5)  Maximum  temperatures  are  shown  to  coincide  with  the  height  of  maximum  leaf  density  (and  minimum  wind 
speeds); 

6)  A  local  minimum  in  temperature  occurs  through  the  trunk  space  below  the  live  branches; 

7)  Under  cloudy  conditions  (with  and  without  rainfall),  mostly  isothermal  conditions  have  been  observed  within 
a  pine  forest  canopy  throughout  the  day; 

8)  At  night,  the  forest  canopy  loses  heat  and  is  cooled  further  (from  the  upper  crown  area  downward)  by 
subsidence  so  that  isothermal  conditions  (or  sometimes  weak  negative  lapse  rate  conditions)  prevail; 

9)  A  second  inversion  is  likely  to  form  close  to  the  ground  at  night  due  to  surface  (radiative)  cooling  of  the  soil; 
and 

10)  Intermittent  wind  gusts  through  openings  at  the  top  of  the  forest  canopy  often  bring  downdrafts  of  cooler  and 
drier  air  from  above  into  the  canopy  and  subcanopy  layers. 
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Table  A -9.  Observed  short-range  acoustic  attenuation  in  forests. 


Parameters 

Reference 

Attenuation  of  sound  waves  in  forests  is  said  to  involve  three  main  phenomena; 

1)  Interference  between  direct  and  ground  reflected  acoustic  waves; 

2)  Scattering  by  tree  trunks  and  branches,  the  ground,  and  turbulence;  and 

3)  Absorption  by  the  trees,  leaves,  branches,  the  ground,  and  the  air 

Ground  impedance,  which  affects  sound-wave  reflection  and  absorption,  is  a 
function  of  flow  resistivity  and  surface  porosity. 

=  20-80  kPa  s  is  flow  resistivity  for  forest  stands; 

=  150-300  kPa  s  mr^  is  flow  resistivity  for  grass; 

=  300-800  kPa  s  mr^  is  flow  resistivity  soil  covered  surfaces. 

Erying  (1946); 

Embleton  (1963); 

Aylor(1972); 

Burns  (1979); 

Embleton  et  al.  (1983); 

Ericke  (1984); 

Attenborough  (1985); 

Price  et  al.  (1988); 

Huismann  and  Attenborough  (1991); 
Sakai  et  al.  (1998,  2001). 

=  0.825  is  the  porosity  for  the  forest  floor,  which  consists  of  humus,  leaves 
(needles),  and  other  biological  litter; 

=  0.675  is  the  porosity  for  grass  fields;  and 
=  0.575  is  the  porosity  for  unpacked  soils. 

•  Ground  impedance  affects  sound  attenuation  mainly  in  the  mid 
250-500  Hz  frequency  range. 

•  At  higher  frequencies  (f  >  1-2  kHz),  sound  attenuation  in  forests  is 
mainly  due  to  scattering  by  tree  trunks  and  branches  and  absorption  by 
foliage. 

•  Such  phenomena  may  have  a  greater  effect  on  sound  attenuation  at 
higher  frequencies  than  refraction  effects  due  to  temperature  inversions 
and  wind  speeds  in  the  forest  canopy. 
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